Evidence for radial flow of thermal dileptons in high-energy nuclear collisions 
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The NA60 experiment at the CERN SPS has studied low-mass dimuon production in 158 AGeV 
In-In collisions. An excess of pairs above the known meson decays has been reported before. We 
now present precision results on the associated transverse momentum spectra. The slope parameter 
Teff extracted from the spectra rises with dimuon mass up to the p, followed by a sudden decline 
above. While the initial rise is consistent with the expectations for radial flow of a hadronic decay 
source, the decline signals a transition to an emission source with much smaller flow. This may well 
represent the flrst direct evidence for thermal radiation of partonic origin in nuclear collisions. 
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Among the observables used for the diagnostics of the 
hot and dense firebaU formed in high-energy nuclear col- 
lisions, lepton pairs are particularly attractive. In con- 
trast to hadrons, they directly probe the entire space- 
time evolution of the fireball and freely escape, undis- 
turbed by final-state interactions. In the invariant mass 
region <1 GeV, thermal dilepton production is largely 
mediated by the broad vector meson p(770). Due to its 
strong coupling to the tttt channel and its short lifetime 
of only 1.3 fm/c, its in- medium properties like mass and 
width have long been considered as the prime signature 
for the restoration of chiral symmetry, associated with 
the QCD phase transition from hadronic to partonic mat- 
ter [ij, y, [Sj . In the mass region > 1 GeV, thermal dilep- 
tons may be produced in either the early partonic or the 
late hadronic phase of the fireball, based here on hadronic 
processes other than tttt annihilation. 

In contrast to real photons, virtual photons decaying 
into lepton pairs can be characterized by two variables, 
mass M and transverse momentum px. Historically, the 
interest has largely focused on mass, including most re- 
cently the first measurement of the space-time averaged 
p spectral function in nuclear collisions [4]. Transverse 
momentum, on the other hand, contains not only contri- 
butions from the spectral function, but encodes the key 
properties of the expanding fireball, temperature and in 
particular transverse (radial) flow. In the description of 
hadron px spectra, the study of collective flow has con- 



tributed decisively to the present understanding of the 
fireball dynamics in nuclear collisions [5|, |6|. However, 
while hadrons always receive the full asymptotic flow 
reached at the moment of decoupling from the flowing 
medium, lepton pairs are continuously emitted during 
the evolution, sensing small flow and high temperature 
at early times, and increasingly larger flow and smaller 
temperatures at later times. The resulting space-time 
folding over the temperature-flow history can be turned 
into a diagnostic tool: the measurement of px spectra 
of lepton pairs potentially offers access to their emission 
region and may thereby differentiate between a hadronic 
and a partonic nature of the emitting source |7|] . 

Experimentally, dilepton px spectra associated with 
direct radiation in the mass region <1 GeV were previ- 
ously only investigated by the CERES/NA45 experiment 
at the CERN SPS (i|, but with low statistics. In this Let- 
ter, we present the first precise data on such spectra, ob- 
tained by the NA60 experiment at the CERN SPS for 158 
AGeV In-In collisions. Mass spectra in different px win- 
dows [9], uncorrected for acceptance, and preliminary re- 
sults on acceptance-corrected px spectra were presented 
before ^Em- 



Details of the NA60 apparatus are contained in |12l.ll3l|. 
The different analysis steps follow the same sequence as 
described in our previous Letter [4], including the crit- 
ical assessment of the combinatorial background from 
TT and K decays by a mixed-event technique [l4J. The 
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FIG. 1: Isolation of an excess above the electromagnetic de- 
cays of neutral mesons (see text). Total data (open circles), 
individual cocktail sources (solid), difference data (thick tri- 
angles), sum of cocktail sources and difference data (dashed). 
Open charm still contained. 

centrality-integrated net mass spectrum of the final data 
sample used here is shown in Fig. [1] It contains about 
430 000 dimuons in the mass range < 1.4 GeV. The nar- 
row vector mesons lo and 4> are completely resolved; the 
mass resolution at the to is 20 MeV. Fig. [T] also contains 
the expected contributions from the electromagnetic de- 
cays of neutral mesons, i.e. the 2-body decays of the ry, lo 
and (j) resonances and the Dalitz decays of the 77, ry and 
uj. The peripheral data can be described by the sum of 
these "cocktail" contributions plus the p 0, Q . This is 
not possible for the total data as plotted in Fig. [3 due to 
the existence of a strong excess of pairs. To isolate this 
excess with a priori unknown characteristics without any 
fits, the cocktail of the decay sources is subtracted from 
the data using local criteria which are solely based on 
the measured mass distribution itself. The p is not sub- 
tracted. The procedure is illustrated in Fig. [1] Each 
centrality window is treated separately. The space of pr 
is subdivided in the range 0<pt<2 GeV into 10 bins of 
equal width, and each bin is also treated separately. The 
yields of the narrow vector mesons w and (j) are fixed so 
as to get, after subtraction, a smooth underlying contin- 
uum, leading to an accuracy of about 2% for the 4> and 
3-4% for the w; at very low pr, even an error of 20% on 
the UJ would still translate to an error of only 5% of the 
excess in the mass region 0.6<M<0.9 GeV. The yield of 
the rj relative to the u) and ^, relevant only for masses 
<0.4 GeV, is fixed from the data at pt>1.0 GeV, profit- 
ing from the very high sensitivity of the spectral shape of 
the Dalitz decay to any underlying admixture from other 
sources. The yield at lower px, subject to lower statis- 
tics due to acceptance, is determined using the NA60 
hadron-decay generator GENESIS |1^] tuned to the cock- 



tail data. The estimated uncertainty on the 77 subtraction 
is ^3%, transforming to a systematic error of ~20% on 
the excess in the mass region 0.2<M<0.4 GeV. Having 
fixed the three main sources, the rj two-body and w Dalitz 
decays are then bound as well; the ratio rj /rj is assumed 
to be 0.12 [15|. Open charm, measured to be 0.30±0.06 
of the total yield in the mass interval 1.2<M<1.4 GeV 
by NA60 [l6|, is subtracted throughout (but not yet in 
Fig. [J), with the spectral shape in M and px as described 
by Pythia [la| ; the resulting uncertainty of the excess is 
8% in this region and <1% everywhere else. After sub- 
traction of the meson decays and charm, the remaining 
sample contains ^150000 dimuons. The subtracted data 
for the J], to and (f) themselves are subject to the same 
further steps as the excess data and are used later for 
comparison. 

In the last step, the data are corrected for the ac- 
ceptance of the NA60 apparatus and for the centrality- 
dependent pair reconstruction efficiencies. The accep- 
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FIG. 2: Transverse momentum spectra of the rho-like mass 
region for three centrality windows, arbitrarily normalized. 
The errors are purely statistical, for systematic errors see text. 

tance shows strong variations with mass and px [Sl, Il0l | , 
but is understood on the level of <10%, mainly based on 
a detailed study of the peripheral data for the particle 
ratios rj/uj and (p/co [J, ISJ . In principle, the acceptance 
correction requires a 3-dimensional grid in M-pT-y space. 
To avoid large statistical errors in low-acceptance bins, 
the correction is performed in 2-dimensional M-px space 
(with 0.1 GeV bins in mass and 0.2 GeV bins in pr), us- 
ing the measured rapidity distribution of the excess as an 
input [lOj . The latter is determined with an acceptance 
correction in y found, in an iterative way, from MC simu- 
lations matched to the data in M and pr- Separately for 
each centrality window, an overlay MC method is used to 
include the effects of pair reconstruction efficiencies. The 
resulting values vary from mostly >0.9 down to about 0.7 
at low mass and low transverse momentum for the high- 
est centrality window. All sources are simulated with a 
uniform distribution in cos^cSj where 9c s is the polar 
angle of the muons in the CoUins-Soper frame, consis- 
tent with (yet unpublished) NA60 data for the w, the (p 
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FIG. 3: Upper: Transverse mass spectra of the excess for 
four mass windows summed over centralities (excluding the 
peripheral bin), in comparison to the </>. Lower: Transverse 
mass spectrum of the rho-like mass region (without uj) for 
very peripheral events, in comparison to the ui. The errors 
are purely statistical, for systematic errors see text. 

and the excess. 

Results on acceptance-corrected pr-spectra for the 
mass window 0.6<A/<0.9 GeV and the three upper cen- 
tralities are shown in Fig. [2j e quiv alent data for other 
mass windows are contained in [10|. Systematic errors, 
not contained in Fig. [21 mainly arise from the uncertain- 
ties of the combinatorial-background and fake-matches 
subtraction and range from 10% to 30% for semiperiph- 
eral up to central collisions at low pt, decreasing rapidly 
to a level of only a few % at higher px- All other uncer- 
tainties are discussed below, in connection with Fig. ID 
The data of the three centrality windows agree within 
their errors; this also holds for the other mass win- 
dows TOJ. 

Fig. [3] (upper) displays the centrality-integrated data 
vs. transverse mass itit, where tot = {Pt + M^)-"^", 
for four mass windows; the (j) is included for comparison. 
The systematic errors due to the background subtraction 
are, at low rriT, 30%, 25%, 15% and 15% in the four win- 
dows, respectively, falling again rapidly to a level of a 
few % at high m^; for the (f>, they are smaller by a factor 
of 5. At very low niT, a steepening is observed in all 
four mass windows, reminiscent of pion spectra and op- 
posite to the expectation for radial flow at masses above 
the pion mass. This increased rise cannot be due to un- 
subtracted background which increases by a factor of 3 
vs. centrality, while the data are independent of it (see 



above). It cannot be due to ill-understood acceptance 
either, since the </>, placed just between the two upper 
mass windows, does not show it. While some rise still 
persists in the peripheral bin, it finally disappears for 
very peripheral coUisions with A<dNch/dri<10 as shown 
in Fig. [3] (lower); the uj does not show it either. The lines 
in Fig. [3] are fits with the function l/rriT dN/drriT oc 
exp{—mT /Tas) , where the effective temperature param- 
eter Tcff characterizes the slope of the distributions. For 
the excess data in Fig. [3] (upper) , the fits are restricted to 
the range 0.4<pT<l-8 GeV (roughly 0.1<(mT-M)<l. 2 
GeV), to exclude the increased rise at low tti^; for all 
other spectra, the fits start at zero. 

The extracted values of Teg vs. pair mass are sum- 
marized in Fig. m (open squares), supplemented by a set 
of further fit values from narrow slices in M (closed tri- 
angles). Preliminary NA60 data from an independent 
analysis [l6| of the intermediate mass region ("IMR") 
1.16<M <2.56 GeV, corrected for the contribution from 
Drell-Yan pairs, are shown for comparison. Finally, the 
hadron data obtained as a by-product of the cocktail sub- 
traction procedure are also included; the value for the 77 
has been obtained by tuning the GENESIS code [15[ to 



the rj Dalitz decay and then referring back to the required 
Teff of the mother. The errors shown for the low-mass 
data ("LMR") are purely statistical. Systematic errors 
only enter to the extent that the slopes of the pt spectra 
are affected, not their absolute level. The errors due to 
the background subtraction are about the same as the 
statistical errors of the fine-bin data (~7 MeV). The er- 
rors associated with the subtraction of the decay sources, 
though significant for the yields (see above), lead to errors 
of only 1-4 MeV for Tcff (dependent on the source), due 
to both the mostly local nature of the subtraction and the 
closeness of Teff for dimuons and hadrons (see Fig.|4|). All 
other error sources considered - relative acceptance, the 
sensitivity to the input y-distribution, cuts vs. no cuts 
in y (none are used), cuts vs. no cuts in cosOcs (none 
are used) - also lead to differences considerably smaller 
than the statistical errors. A correction for Drell-Yan, 
using an extrapolation down to M<1 GeV [TTJ, would 
systematically lower the values by 5-10 MeV, depending 
on mass. 

The results displayed in Fig. [4] can be summarized and 
interpreted as follows. The slope parameter Teg rises 
nearly linearly with mass up to about 270 MeV at the 
pole position of the p, followed by a sudden decline to 
values of 190-200 MeV for masses >1 GeV. The excess 
yield in the mass region 0.2<M<0.9 GeV has generally 
been attributed to thermal radiation from the fireball, 
dominated by pion annihilation tt+tt^ ^^ P ^ M^/^ i 
and the NA60 data, before acceptance correction, have 
directly been interpreted as the space-time averaged in- 
medium spectral function of the p [4] . Following earlier 
work [3|, they are now nearly quantitatively described 
by the newest theoretical developments [171 llSl |l9| . The 
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FIG. 4: Inverse slope parameter Tas vs. dimuon mass M for 
dNch/dri>30. The open squares correspond to the fit hnes in 
Fig. |3] Open charm is subtracted throughout. For explana- 
tion of the inserted symbols and the errors see text. 

linear rise of T^s with M over the whole region up to the 
p peak is reminiscent of radial flow of a hadronic source. 
While the hadron data show a similar linear rise, their 
absolute values are surprisingly close to the excess val- 
ues, contrary to the expectation for the temperature-flow 
folding expressed in the introduction. The seeming con- 
tradiction can be resolved by comparing the (free) p itself 
rather than the other hadrons with the in-medium emis- 
sion. This p is accessible as the peak on the broad con- 
tinuum (Fig. [T]) , generally interpreted as the freeze-out p 
without in-medium effects 17|, llSl . Il9l | . By disentangling 
the peak from the continuum as done before [9|, llO| , we 
find Toff=300±17 MeV for the peak and 231±7 for the 
underlying continuum in the window 0.6<M<0.9 MeV. 
The high value of the peak, added as a further hadron 
point into Fig. m should then be interpreted as charac- 
teristic for the true freeze-out parameters of the fireball, 
implying the ?/, lu and cf) to freeze out earlier, due to their 
smaller coupling to the pions. By modeling a p with this 
temperature, its contribution can be subtracted from the 
total for each of the finer binned data points in Fig. IH 
lowering Tes by 4—20 MeV depending on the closeness to 
the pole, and shifting the maximum of the resulting "in- 
medium" values up to the mass bin just below 1 GeV. In 
any case, the large gap in T^s between the vacuum p and 
the excess points (corrected or not for the subtraction) 
restores the expected difference between a freely emitted 
hadron and its in-medium decay part, making the ob- 
served linear rise of Teg with M now consistent with the 
expectations for radial flow of a hadronic source (here 
TTTT — > p) decaying into lepton pair s. Theoretical model- 
ing of our results is underway [l7|, ll8J, |l9| , but does not 
yet describe the data in a satisfactory way. 

It is interesting to note that the large gap of >50 



MeV in Teff between the vacuum p and the lu (same 
mass) decreases towards the peripheral window, but 
only closes for the lowest peripheral "pp-like" selection 
4:<dNch/dr]<10 shown in Fig. [3] (lower), with roff=198±6 
MeV for the p and 201±4 MeV for the w. This implies 
that both the "hot p" and the low-m^ rise discussed be- 
fore are intimately connected to pions, disappearing to- 
gether as the TTTT contribution to p production vanishes 
(with only the cocktail p left). 

The sudden decline of Toff at masses >1 GeV is the 
other most remarkable feature of the present data. Ex- 
trapolating the lower-mass trend to beyond 1 GeV, a 
jump by about 50 MeV down to a low-flow situation is 
extremely hard to reconcile with emission sources which 
continue to be of dominantly hadronic origin in this re- 
gion. // the rise is due to flow, the sudden loss of flow 
is most naturally explained as a transition to a quali- 
tatively different source, implying dominantly early, i.e. 
partonic processes like qq — > p.'^ p~ for which flow has not 
yet built up [l8|. While still controversial [17|], this may 
well represent the first direct evidence for thermal radi- 
ation of partonic origin, breaking parton-hadron duality 
for the yield description in the mass domain. 

In conclusion, we have found strong evidence for radial 
flow in the region of thermal dilepton emission which has 
previously been associated with the p spectral function. 
The transition to a low-flow region above may signal a 
transition from a hadronic to a partonic source. 
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